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Abstract 

I propose event-plane correlations as a test of collectivity in small colliding systems: In d-Au 
and 3 He-Au collisions, I predict a strong anti-correlation between elliptic and triangular flow, 
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I. INTRODUCTION 


In high energy nucleus-nucleus collisions, collectivity is best manifested in the measured 
harmonic flow [1, 2], Defined as the Fourier harmonics of a single-particle spectrum, 

-J— = + Vne~ m<t>p + complex conj.], (1.1) 

d(b v 2 tt 

y n =1 

harmonic flow V n = v n exp(m'h„) characterises final state anisotropy in the momentum 
space. Recently, analyses of harmonic flow were extended to small colliding systems, re¬ 
garding the striking observation of a ridge structure in long-range (in rapidity) two-particle 
correlations in high multiplicity events [3] . At RHIC energy, substantial elliptic flow v 2 and 
triangular flow v 3 were extracted in d-Au and 3 He-Au collisions [4, 5]. At the LHC, in 
^/s NN = 5.07 TeV p-Pb collisions, it was also shown that harmonic flow can be as large as 
that measured in Pb-Pb collisions with = 2.76 TeV [6-8], provided that multiplicity 
yields in these events are comparable [9] . The observation of harmonic flow implies a possible 
collective expansion stage during the collisions of p-Pb, d-Au and 3 He-Au, which is further 
supported by various theoretical simulations. With appropriate descriptions of the initial 
density profile, reasonable predictions on harmonic flow have been achieved for p-Pb [10- 
12], d-Au [12, 13] and 3 He-Au [13, 14] colliding systems with viscous hydrodynamics. In 
a transport approach, it is also found that the prediction of v 2 and v 3 requires sufficient 
interactions in the late stage [15]. 

However, the statement of collectivity in small colliding systems is still under debate. In 
particular, it is challenged by the fact that long-range correlations seen in the experimental 
data can originate alternatively from pre-equilibrium physics [16]. A way to disentangle 
competing models, and especially to clarify the issue of collectivity, lies in the measurements 
of flow fluctuations and mixing between harmonics. The idea relies on the knowledge that 
medium collectivity associates harmonic flow with event-by-event fluctuations in the initial 
density profile. The CMS Collaboration has updated its measurements of v 2 fluctuations in 
terms of cumulants v 2 {m} from multi-particle correlations [17, 18]. The observed pattern of 
v 2 {m} agrees with cumulants of initial cllipticity which is purely driven by fluctuations [19, 
20], which coincides with the expectation from collective expansion. 

This work is motivated in a similar manner, but with emphasis on the mixing of harmonic 
flow. Especially, event-by-event fluctuations in the initial density profile dominate the mixing 
between lower harmonics due to a linear eccentricity scaling [21, 22]. This allows one to 
estimate the mixing between V 2 and V 3 for p-Pb, d-Au and 3 He-Au, in terms of the mixing 
between initial cllipticity and triangularity. Analyses in this work are focused on central 
collision events, where the initial state geometries differ dramatically in p-Pb, d-Au and 
3 He-Au, with the background density profile azimuthally symmetric, dumbbell-shaped and 
triangle-shaped respectively. In addition, correlations of harmonics involving a dipolar flow 
V\ are studied as well, based on initial mixing involving dipolar asymmetry. 

In order to have a consistent comparison with experimental measurements, a series of cor¬ 
relation coefficients are specified in Section II. Collision events are simulated by PHOBOS 
Monte Carlo Glauber model, with details of the model simulation described in Section III. 
Correlations of initial anisotropies from simulations are presented in Section III. In Sec¬ 
tion IV, analytical analysis of the correlations is given based on an independent source 
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approach, where initial state of the small colliding systems are modelled by a number of 
independent sources. Effects beyond linear eccentricity scaling are discussed with respect to 
the mixing between V 2 and V 3 in p-Pb collisions in Section V. 


II. CORRELATION COEFFICIENTS 


In this paper, to quantify correlations between V n and V m (i.e., phase mixing between 
and when n 7 ^ m), we define the following correlation coefficient which has a similar 
structure as a Pearson correlation coefficient, 

v „ Mv:(v,i) b ) ,,,, 

nm (iKi 2 “) i/2 <iu.n i/2 ’ K -’ 

where a and b are integer numbers, satisfying ma — nb — 0 due to a rotational symmetry. 
In Eq. (2.1) and throughout this paper, angular brackets (...) are used to notate average 
over events. By definition one has 0 < V^m < 1, with 1 stands for an absolute correlation 
or anti-correlation and 0 corresponds to pure randomness. 

has been studied through the measurements of event-plane correlations in Pb-Pb 
systems by the ATLAS collaboration [ 21 ], using the scalar-product method. In particular, 
the correlation between V 2 and V 3 , 


Re(Vi(V 3 ’) 2 ) 

23 <| V 2 | 6 > V 2 (| V ' 3 | 4 ) 1 / 2 ’ 


( 2 . 2 ) 


was found negligibly weak for most of the centralities, except peripheral collisions. Unlike 
correlations involving higher harmonics, which receive signihcant contributions from non¬ 
linear flow response during medium expansion [ 22 ] , is mostly determined by the mixing 
between initial anisotropies. Anisotropies of initial state are quantified generally by eccen¬ 
tricities. In each single collision event, if one denotes average over transverse plane according 
to the initial density profile as {...}, the Erst three eccentricities are 1 


£, =e 1 e i$1 = - 
£2 =e 2 e i2 ^ = 
£3 =£se i3 * 3 = 


'“FT’ 

2^ r 

{ r 2} 

{ r 3} 


(2.3) 

(2.4) 

(2.5) 


In accordance with Eq. (1.1), complex notations have been applied in Eqs. (2.3), (2.4) 
and (2.5), so that both magnitude and phase of eccentricity are defined simultaneously. 
Approximately a linear eccentricity scaling can be assumed for lower harmonics, based on 
event-by-event hydrodynamic simulations [23], 


Vk K n £ n , Tl ^ 3 . 


( 2 . 6 ) 


1 Re-centering corrections are implied in these definitions, with {re l ^ r } = 0 
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K n in Eq. (2.6) is known as the medium response coefficient, which is determined by the 
property of medium collectivity. It is then obvious that correlations of final state between 
harmonics are identical to the corresponding correlations of initial anisotropies, e.g., 

^ 23 =^ 3 - ( 2 - 7 ) 

Note that the superscript with e’s on the right hand side of Eq. (2.7) implies correlations of 
initial anisotropies. 

Mixing between S 2 and £3 (or equivalently V 2 and V 3 ) in small colliding systems is of 
particular interest for two reasons. First, geometrical fluctuations in the initial state are 
much more pronounced in p-Pb, d-Au and 3 He-Au than those in nucleus-nucleus collisions. 
For instance, the multiplicity in central p-Pb is smaller roughly by a factor of 10 than in 
central Pb-Pb collisions, one expects initial state fluctuations to be larger by a factor of 
VTO. Second, background geometries in d-Au and 3 He-Au are generically deformed. For 
events in the most central collision bins of p-Pb, d-Au and 3 He-Au, geometry of the system 
is dominated by the configuration of proton, deuteron and 3 He respectively. Therefore, 
initial state background geometry of p-Pb is azimuthally symmetric. But for d-Au and 
3 He-Au, a dumbbell-shaped initial density and an intrinsic triangle-shaped initial density 
are expected respectively. Accordingly, non-trivial geometries lead to an intrinsic initial 
cllipticity £ 2 in d-Au and an intrinsic initial triangularity £3 in 3 He-Au. Both fluctuations and 
background geometry will be shown essential for the generation of anisotropy correlations, 
later in Section III and Section IV. Similarly, one also expects fluctuations and geometric 
shape induce correlations involving a dipolar asymmetry, namely correlations of harmonics 
involving V±. In this paper the following types of correlations are studied as well, 

Vr 2 , V\l and VIH- ( 2 . 8 ) 

Note that is defined as 

= Re jffigj 

123 (lUUlWN 2 ) 1 / 2 ’ { ’ 

which is normalized differently from the three-plane event-plane correlations introduced by 
the ATLAS collaboration [ 21 ] with the scalar-product method. 


III. MONTE CARLO GLAUBER SIMULATIONS 


In this work, we simulate initial state of heavy-ion collisions using the PHOBOS Monte 
Carlo Glauber model [24, 25]. Additional modifications are taken into account with respect 
to small colliding systems. For instance, the wounding profile of two colliding nucleons is 
taken to be a Gaussian instead of a step function [26]. Regarding heavy-ion collisions, PHO¬ 
BOS Glauber model captures the dominant geometric structure and effects of fluctuations 
of nucleons inside a colliding nucleus. On an event-by-event basis, fluctuations of nucleon 
positioning inside nucleus are introduced randomly in model simulations, resulting in a dis¬ 
crete distribution of number of participants N part and collisions N co u. A Gaussian smearing 
centered at the location of each participant and collision is applied, 


1 -IUJW 
27TCT 2 


(3.1) 
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FIG. 1. (Color online) Initial state mixing between £2 and £3 as a function of number of par¬ 
ticipants. Results in panel (a) are from PHOBOS Monte Carlo Glauber model with a constant 
weighting of C s in Eq. (3.2). Dashed line is analytical expectation from the independent-source 
model. Results with an extra NBD weighting is shown in Panel (b). Results from smearing with 
a = 0.4 frn are plotted with symbols, while shaded bands correspond to the case with a = 0.8 fm. 


The size of smearing is taken to be cr = 0.4 fm and 0.8 fm in this work for analysis. 
Initial entropy density is thus obtained by assuming a proper weight associated with each 
participant and collision, 


s(x,y,r 0 ) = C s 


1 



N P art(x , y ) + aN cM (x, y) 


(3.2) 


Instead of a constant weighting of C S) a negative binomial distribution (NBD) weighting is 
known necessary to reproduce the measured multiplicity distributions [26] in small colliding 
systems, which introduces extra fluctuations. Especially, probability of generating high 
multiplicity events is enhanced with a negative binomial weighting. In this work, following 
Refs. [12, 26] and taking the same set of parameters for p-Pb at the LHC, and d-Au and 
3 He-Au at RHIC, ten million events are generated. Events are then classified with respect 
to number of participants N part , which has a correspondence to the centrality classification 
used in experiments. For simulations with ten million events, the centrality of bin class with 
reliable statistics can be as small as 0.01%, which corresponds to the right most points in 
Fig. 1 and Fig. 2. 

Anti-correlations between S 2 and £3 in p-Pb, d-Au and 3 He-Au are discovered via model 
simulations, with results shown in Fig. 1. Note that a minus sign is implied. Panel (a) of 
Fig. 1 is obtained based on PHOBOS Glauber model, without additional NBD weighting 
for the entropy production. Albeit less realistic, the case with a constant weighting is 
closer to the model of independent sources [27? , 28], especially for the p-Pb system the 
correlation "Pff is found quantitatively consistent with predictions from independent sources 
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FIG. 2. (Color online) Initial state mixing involving £\ as a function of number of participants, 
represented by (a) 'Pff, (b) -"Pff, and (c) -Vf^- All the results are obtained with an extra NBD 
weighting for the entropy production. Results with a = 0.4 fm are plotted with symbols, while 
shaded bands correspond to cr = 0.8 fm. 


(see discussions around Eq. (4.3) in Section IV ). 2 Simulations with NBD weighting which 
correspond to experiments carried out at RHIC and the LHC are presented in Fig. 1 (b). 
Errors of these results are estimated based on variance. Effect of smearing is investigated 
with a = 0.4 fm (symbols with error bars) and cr = 0.8 fm (shaded bands). Larger smearing 
size leads to stronger correlations between £2 and £ 3 , which is more significant for d-Au and 
3 He-Au systems. 

In all three colliding systems, the correlation between £2 and £3 decreases towards central 
collisions. This decrease is associated with the decrease of fluctuations in the initial state. 
When extra fluctuations are introduced from a NBD weighting of entropy production, seen 
as changes from Fig. 1 (a) to Fig. 1 (b), correlations are enhanced as well. Especially for 
p-Pb, the correlation of the most central events reaches 10%, as shown in Fig. 1 (b). In 
addition, from Fig. 1 there is an apparent hierarchy relation of correlations among the three 
systems, 

ra P -Pb)| < | 7 >f 3 £ ( 3 He-Au)| < |Rg(d-Au)| (3.3) 

which reflects the differences in the background geometry of each system. One may naively 
understand this relation by noticing that £2 has one power higher than £3 in the original 
definition of Pff, which makes a backgound cllipticitiy more important than a background 
triangularity. More detailed discussion of this hierarchy relation are given in Section IV. 

Fig. 2 presents initial correlations involving £±, as predictions for the flow correlations 
involving V). Again, to have all results on the same scale a minus sign is applied when 
necessary. Irrespective of the detailed background geometry, in all the three systems a 
negative and a negative are obtained which indicate anti-correlations between £\ 
and £ 2 , and among £\, £2 and £3 respectively, while £\ and £3 are found strongly correlated. 

2 Note that in the original PHOBOS Glauber model[24, 25], a repulsive correlation is applied in the 
simulations, so it is not strictly equivalent to a model of independent sources. 
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FIG. 3. (Color online) Alignment of anisotropies due to fluctuations on top of smooth background 
(solid blue shape) with a net ellipticity £ 2 (a) and a net triangularity £3 (b). Green dashed shapes 
correspond to anisotropies generated by one addition hot spot in the tip area of the background. 
Dipolar asymmetry is illustrated as green arrows. 


Although these correlations are stronger than 'Pff, they exhibit very similar dependence on 
the background geometry of the system and fluctuations, i.e., similar hierarchy relations 
regarding the three systems. More detailed discussions are postponed to Section IV. It is 
also worth mentioning that for the most central events in d-Au systems, correlations from 
simulations are compatible with those expected in peripheral nucleus-nucleus collisions [29— 
31] , where the effects of fluctuations are comparable. 

Measuring event-plane correlations in small collision systems is a challenging analysis. In 
Appendix A, we evaluate the order of magnitude of the statistical error. We find that in 
order to detect with a few percent accuracy in p-Pb central collisions at the LHC, one 
typically needs a million events per centrality bin. 


IV. UNDERSTANDING OF INITIAL STATE CORRELATIONS 

The origin of anisotropy correlations in the initial state of heavy-ion collisions, with 
respect to the effects of fluctuations and a background geometry, can be understood qualita¬ 
tively as follows. In heavy-ion collisions, a bumpy initial state density profile can be equiva¬ 
lently realized by randomly throwing additional spots on top of a smooth background. For 
a smooth background with a net ellipticity or triangularity, the area of tip is always larger 
than that of side, which is purely a consequence of geometric effect. Therefore, one expects 
the additional spots to appear more probably in the area of tip rather than side. Fig. 3 
describes anisotropies induced by one additional hot spot in the tip area of an elliptic back¬ 
ground (a), and a triangular background (b). In both cases, one can check that the induced 
anisotropies are aligned with the background shape in a specified way: $1 = tt, $ 2 = 7t/ 2, 
and $3 = 7 t/ 3 , which generates perfect anti-correlations between £ 2 and £ 3 , S\ and £ 2 , and 
among £ 1 , £ 2 and £ 3 , and perfect correlation between S\ and £ 3 . For a more sophisticated 
density profile which fluctuates from event to event, correlations are weaken due to contri¬ 
butions from configurations with excessive density around sides of the background shape, 
but the sign of correlation should not be changed. 
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To a quantitative level, the effects of fluctuations and background geometry on correla¬ 
tions can be studied in the independent-source model [27], in which an initial state density 
profile is modelled by N point-like independent sources on top of a specified background. 
Consequently, fluctuations of any quantity / is characterized as 


8 f = {/} - (/), 


(4.1) 


where {...} and (...) are used in the similar way as in Eqs. (2.3) and (2.1), indicating 
average in the transverse plane for one single event and average over events respectively. 
Independency of sources allows one to write event averaged quantities in powers of 1/IV, 
i.e., fluctuations. For instance, the two-point function which quantifies event average of 
quadratic order of fluctuations, is of the order of 1/IV, 

/r r \ _ (fa) - ( f)(a > u 0 x 

(°f°g) ~ jy ■ ( 4 - 2 ) 

Event average of higher order of fluctuations leads to higher order dependence on 1/N. More 
details of the model can be found in [27] and in Appendix B. For p-Pb collisions, background 
is azimuthally symmetric and thus £ 2 = £3 = 0. One accordingly finds that the dominant 
contribution is purely driven by fluctuations, and of the order of 1/1V 3 / 2 , 


Pg(p-Pb) 


J_ (r 12 ) + 6{r fi ) 2 — 18(r 4 )(r 8 ) 5/2 

iV s / 2 Vl2{r^)s/ 2 {r«) 1 ' ’ 


(4.3) 


For a Gaussian smearing, the geometric factor in Eq. (4.3) is known explicitly (« 1.22) and 
thus an analytical description can be obtained. As shown as the blue dashed line in Fig. 1 
(a), the prediction from independent sources naturally captures the correlation pattern if one 
assumes that the number of independent sources scales like N part ~ 0.55IV. Note also that 
in Eq. (4.3) the correlation is insensitive to a Gaussian smearing, which is indeed confirmed 
in Fig. 1 that smearing a has little influence on the simulated correlation 'Pff(p-Pb). 

In principle, the generic correlations induced by the nuclear configurations of deuteron 
and 3 He prevent applying the independent-source model to the collisions of d-Au and 3 He- 
Au. Nevertheless, the effect of a non-trivial background geometry is worth analyzing, at least 
qualitatively. For 3 He-Au collisions, the conditions £2 = 0 and £3 7 ^ 0 allow terms of lower 
order in 1/IV contribute. While for d-Au systems, a non-zero background ellipticity picks 
up even lower order terms in the expansion with respect to 1/IV. It can also be understood 
directly from the definition of "Pff ~ (^K^g) 2 ), that £2 has one power higher than £ 3 . In 
total, counting the extra suppression due to fluctuations, one finds 


|P| 3 £ ( 3 He-Au)| ~ O 



1^23(d-Au)| ~ 0(1), 


(4.4) 


which is consistent with Eq. (3.3). It is also interesting to notice that the hierarchy relation 
in Eq. (3.3) disappears when N part is sufficiently small, as depicted in Fig. 1, where only one 
of the constituent nucleons in deuteron and 3 He participates the collision. 

For correlations involving a dipolar asymmetry, an ideal assumption would be taking £\ 
as a result of fluctuations in all the three small systems, then the similar argument leads to 

|GL 

Wn 



|Pffi(d-Au)| ~ |T>S( 3 He-Au)| ~ 0(1), |^S(p-Pb)| ~ O 


(4.5) 






due to the fact that £2 and £3 are of the same power in the definition of 'Pfff. And indeed 
Eq. (4.5) agrees with the results depicted in Fig. 2 (c). However, when the role of £\ becomes 
more important, as in "Pff and Pff which contain £f and £f respectively, one has to take into 
account the effect of an intrinsic £\ in both d-Au and 3 He-An 3 . Because a larger intrinsic 
£\ is expected with repect to a dumbbell-shape than a triangular shape, the correlations 
and are stronger in d-Au than in 3 He-Au. 


V. IN PROTON-LEAD BEYOND LINEAR ECCENTRICITY SCALING 


One must go beyond the approximate linear eccentricity scaling to have a more realistic es¬ 
timate for the correlations of harmonics in experiments. For the central p-Pb collisions, non¬ 
linear flow response can still be ignored since background geometry is azimuthally symmetric, 
which case is similar to the ultra-central Pb-Pb. The most significant corrections come from 
fluctuations in hydro response. While detailed knowledge of event-by-event fluctuations of 
flow response requires a systematic analysis based on, e.g., event-by-event hydrodynamic 
simulations, it can be quantified by introducing an extra noise term in Eq. (2.6) [32, 33] 

V n = K n £ n + X n , (5-1) 


where the noise X n is complex. I further model X n as a Gaussian noise, which is uncorrelated 
with the initial eccentricity. This noise breaks linear eccentricity scaling [23]. The magnitude 
of the effect can be quantified by the Pearson correlation between the flow and the initial 
anisotropy, 


(Vn£S 

<|l4| 2 }l/2(|£„|2)l/2 


1 — A n + 0{ A;), 


(5.2) 


where 


1 (|AA| 2 ) 

2^<|£nir 


(5.3) 


is a small and positive quantity, characterizing the relative magnitude between fluctuations 
and average flow response. For the correlation between V 2 and V3, substituting Eq. (5.1) 
into Eq. (2.2) results in 


V vv 

1 23 


j>ee 
' 23 



9 (|£| l )(|£ 2 | 2 ) 

<|£ 2 | 6 > 


As 


4<|£ 3 | 2 ) 2 

<|£ 3 | 4 > 


A 3 + 0( A 2 ) 


(5.4) 


The factors before A 2 and A 3 in Eq. (5.4) are determined by fluctuations of £ 2 and £3 
respectively. In p-Pb collisions, it was proposed that fluctuations of initial state eccentricities 
follow the so-called power distribution [19], from which the ratios between moments of initial 
state eccentricity are known explicitly. In central collision events, approximately one has 
9<|£ 4 |)<|£ 2 | 2 >/<|£ 2 | 6 > * 3 and 4(|£ 3 | 2 ) 2 /(|^| 4 ) « 2. 

An estimate of A 2 and A 3 in p-Pb collisions at the LHC can be made in terms of the 
breaking of two-particle correlations. Fluctuations in hydro response result in a breaking in 
the two-particle correlation function, which is characterised by the ratio [34], 

r 


1 (PtjPt) ~ 


An A {JPti Pt ) 




(5.5) 


3 Note that there is no intrinsic £\ in p-Pb. 
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Taking into account noise in hydro response in Eq. (5.1), one obtains 

rniPTiPr) = 1 - K(Pt) ~ K(Pt) + °(^ 2 )- (5.6) 

For the central collision events in the available pp range, ^(p^j-Pt) from event-by-event 
hydrodynamic simulations can reach as low as 97% [11], depending on smearing size and 
shear viscosity over entropy ratio p/s used in simulations. Consequently, one would expect 
that the integrated value of A 2 to be smaller than 3 %. 4 Similarly, the upper bound of A 3 
can be taken as 5%. 

In total, for the central events of p-Pb, correlation between V 2 and V 3 is expected with a 
value greater than 8 % if late stage collective expansion indeed dominates the evolution. 


VI. CONCLUSIONS 

In summary, I have studied correlations between harmonics in the small colliding systems: 
p-Pb, d-Au and 3 He-Au. Emphasis is laid on the mixing between V 2 and V 3 , in terms of 
mixing between cllipticity £2 and triangularity £ 3 , by assuming linear eccentricity scaling. 
Significant anti-correlations are found between £0 and £3 in p-Pb systems from Monte Carlo 
Glauber simulations, with |'Pff| of the order of 10% in the central collision events. When 
the background geometry exhibits an intrinsic ellipticity, such as central d-Au collisions, or 
an intrinsic triangularity, such as central 3 He-Au collisions, the correlation is even enhanced. 
Monte Carlo Glauber simulations present stronger anti-correlation of "Pff in d-Au than that 
in 3 He-Au. Similar, but stronger correlation patterns are found when the dipolar asymmetry 
£\ is involved. The physical origin of initial state correlations is analytically studied in the 
independent-source model, from which fluctuations and intrinsic shape are recognized as the 
dominant effects. 

Mixing between harmonic flow is basically considered equal to the corresponding mixing 
of initial anisotropies, due to an approximate linear eccentricity scaling. Effects beyond 
linear eccentricity scaling are discussed by including noises in the linear flow response, from 
which a more realistic estimate of V 23 is made for the central p-Pb collisions. Regarding 
measurements carried out by the CMS collaboration, statistics in experiment is discussed in 
Appendix A. Mixing between V 2 and V 3 is thus proposed as an accessible probe to detect 
the medium collectivity in experiments, if evolution in p-Pb collisions is dominated by the 
medium collective expansion. In d-Au and 3 He-Au, especially in the most central events, 
medium collectivity is supposed to lead to sizable as well, although its patterns may 
be more involved due to other effects like non-linearities in the flow generation. However, 
geometries of p-Pb, d-Au and 3 He-Au systems are expected to be revealed by comparing V 23 
among these colliding systems. I 11 a similar manner, correlations involving V\ are proposed 
as probes as well, which are generated from initial correlations involving £\. 


4 Actually, to obtain the value in Eq. (5.3) the integration with respect to pr should be done separately for 
numerator (f dpr(\X n (pT)\ 2 )dN/dpr) and denominator (f dpT^K^p^H^l^dN/dpr) , which leads to an 
even smaller ratio than f dpT^ 2 {pr)dN / dpr ■ 
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Appendix A: Statistical errors of measurements in p-Pb 

In this appendix, I briefly estimate the statistical error on in the central p-Pb col¬ 
lisions at the LHC. The order of magnitude of the number of events that is needed for the 
observation of V 23 in experiments can be inferred consequently. The measurement of "P 23 
involves a five-particle correlation, 

<E 2 3 (H 3 *) 2 > (Al) 

whose statistical error is determined by the total number of independent 5-plets which can 
be constructed. For a total number of events N eve and multiplicity M in each single event, 
statistical fluctuations with respect to (V^ 3 (V^*) 2 ) are therefore of the order 1/ (N eve M b ) 1 / 2 . 
Statistical errors on V% 3 from measurements can be accordingly estimated as 

YP™ ~ —L= • (A2) 

VN^e \v 2 Vm) \v 3 Vm) 

The rule of thumb of writing Eq. (A2) is that each V n contributes a factor 1/(|V^| \/M). Sim¬ 
ilar analysis can be generalized to the measurements of other types of correlations between 
harmonics. 

Referring to the measurements of p-Pb carried out by the CMS collaboration [6, 18], the 
resolution of v 2 and v 3 measurements for central collisions are approximately (1/^2 \fM) ~ 
1.25 and (1 / v 3 \[M ) rs./ 5 respectively. In total, to observe P 23 with value at least 8%, a 
total number of one million events are roughly needed for each centrality bin. The LHC has 
already recorded an integrated luminosity of 35nb _1 (cf. Ref. [18]) in the p-Pb runs collected 
in 2012 and 2013, corresponding to a total events of the order of 10 10 . Therefore there are 
over one million events in the 0.01% most central bin, and the anti-correlation between V 2 
and V 3 can be detected as a feasible probe for the test of medium collectivity. 


Appendix B: P 23 from independent sources 


We take complex notations for simplicity in the following derivations, so that the trans¬ 
verse coordinate is expressed as z — r exp In the independent-source model, azimuthal 

eccentricities of a fluctuating initial state can be analytically written. Taking into account 
of the re-centering corrections, 


£ 2 —■^2 e 


*2$ 2 


C —c 

£3 ~ £ 3 e 


{{z-S*) 2 } 

{ r 2} 

{{z-m 

{r 3 } 


(Bl) 
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Eq. (Bl) can be expanded with respect to fluctuations order by order, which in turn gives 
rise to an expansion in terms of 1/N. Therefore we have, 


£2 

£3 


[( z 2 ) + 5 z 2 + 

p> 


+ o(£ 3 ), 


[(z 3 ) - 3(z 2 )6 z + S z 3 - Z5 z 5 z i] 3 


( r 3) 


+ 0(V 


(B2) 


Note that ( z 2 ) = — (r 2 )£2 or ( z 3 ) = — {r 3 )£^, in Eq. (B2) vanishes when there is no corre¬ 
sponding intrinsic eccentricity from the background geometry. For simplicity, contributions 
from fluctuations of denominators in Eq. (Bl) are neglected, which does not affect the 
counting of power of fluctuations and is allowed quantitatively for p-Pb. 

In the independent source model, in addition to the two-point function in Eq. (4.2) which 
is proportional to 1/N, multi-point function can be evaluated using Wick’s theorem, and 
found to be suppressed by higher powers of 1/N [35]. Implicitly, one has 


(■5 3 > ~ <U ~ O (- 4 ) . (5 3 ) ~ (5‘) ~ O (A) ,... . 


(B3) 


Pearson coefficients of initial state Eq. (2.1) can therefore be written order by order in powers 
of 1/N. For the p-Pb systems, because all terms proportional to (z 2 ) or (z 3 ) vanish, one has 


V ££ — 

' 23 — 


(f>%5%) + 3{5%5 2 S%) - 6 (5%5A*S„) + 0((5 7 )) 


[<%->> + 0 ((^))] 1/2 [(5%5%) + 0««'))] 


1/2 


iV 3 / 2 


(B4) 


The minus sign in Eq. (B4) comes from convention used in the definition of anisotropies and 
which agrees with the analysis given in Section IV. Again, in evaluating multi-point 
functions, azimuthal symmetry in p-Pb demands that only terms of the type (r 2n ) = ( z n z n ) 
left, which leads to the result in Eq. (4.3). For 3 He-Au, terms containing (z 3 ) contribute. 
Besides, terms of the type (z a z b ) 7 ^ 0 also contribute if | a — 6|/3 = 0, 1, etc., which leads to 


<.° 3 ) 2 (fe) + 3(s 3 ) 2 (^g) + + om) 

{(z^(z*y + 0«<5 2 »p 2 [(51,5%) + 0«<5«»] 1/2 

(z s ) 2 (z e )+6(z% 2 (z 3 ) 2 t 1 \ 

VW(z 3 )(z 3 ) (r ‘) 3 / 2 + U V N s/2 ) ' 


For d-Au collisions, (z 2 ) 7 ^ 0 and accordingly all terms of the type ( z a z b ) 7 ^ 0 if \a — b\ is a 
even integer. One finds that the lowest order term is independent of 1/iV, 
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